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Determination of as from Event Shapes and Power Corrections 

Daniel Wicke'* 

'^Fachbereich Physik, Bergische Universitat GH, Gaufistr. 20, 42097 Wuppertal; wicke@cem.cli 

The size of non-perturbative corrections to event shape observables is predicted to fall like a power of the 
inverse centre of mass energy. These power corrections are investigated for different observables from e^e~-anni- 
hilation and compared to the theoretical predictions. Using the latest DELPHI high energy data advantages of 
determining as from these predictions are discussed and compared to conventional methods. 



1. Introduction 

The upgrade of LEP to run at energies above 
the Z-peak induced new possibilities in measur- 
ing the running of as- Due to asymptotic free- 
dom, non-perturbative effects become less im- 
portant at higher energies. Different theoretical 
frameworks suggest that for event shapes the non- 
perturbative corrections reduce like the inverse 
power of the centre of mass energy Theo- 
retical predictions for this decrease can be used 
in the measurement of independent of event 
by event hadronization models. 

2. as from Mean Event Shapes using 
Power Corrections 

The analytical power ansatz for non-perturba- 
tive corrections by Dokshitzer and Webber 
can be used to determine as from mean event 
shapes This ansatz provides an additive 

term to the perturbative 0{a'^) QCD prediction 

(/pert) ■ 



{f) = — f f^da = (/pert) + (/pow) 

CTtot J da 
The power correction is given by 
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an infrared matching scale /^/, K — (67/18 — 
7rV6)CA - 5A^//9, fi^ = (33 - 2Nf)/12TT and 
a,f = 'iCf/n. Beside as this formulae contains 
ao as the only free parameter. In order to mea- 
sure as this parameter has to be known. 

To infer ao a combined fit of as and ao to a 
large set of measurements at different energies 
is performed. For Ecm > Mz only DEL- 
PHI measurements are included in the fit. Fig. || 
shows the measured mean values of (1 — T) and 
(M^/E'^jg) as a function of the centre of mass 
energy together with the results of the fit. The 
resulting values of ao are summarized in Tab. 
The extracted ao values are around 0.5 as ex- 
pected in 1^, but they are incompatible with 
each other. As the assumed universality is not 
given to the precision that is accessible from 
data, ao is individually determined for (1 — T) 
and (Af^/i?^ig). The scale error is obtained from 
varying the renormalization scale / — ^'^/E^^ 
from 0.25 to 4 and the infrared matching scale 
from IGeV to 2GeV. 

After fixing ao, the as values corresponding 
to the high energy data points can be calculated 
from Eqs. (0-1). as is calculated for both ob- 
servables individually and then combined with an 



Table 1 

ao as determined from mean event shapes. 



ao is a non-perturbative parameter accounting 
for the contributions to the event shape below 



Observable 


ao ± Stat. ± scale 


X^/ndf 


(1-T) 
(Mi/E^^,^) 


0.531 ±0.012 ±0.003 
0.434 ±0.010 ±0.010 


42/23 
4.0/14 
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Table 2 

as as obtained from means with the Eqs. 

by averaging the (1 — T) and (M^/E'^j^) resuhs. 





oi-s{Eci-n) ± Stat. ± sys. ± scale 


133GeV 
leiGeV 
172GeV 


0.1162 ±0.0070 ±0.0009 ±0.0047 
0.1047 ±0.0069 ±0.0034 ±0.0043 
0.1052 ±0.0082 ±0.0037 ±0.0041 



unweighted average. Its error is propagated from 
the data and combined by assuming full correla- 
tion. An additional scale error was calculated by 
varying / and /i/ in the ranges discussed. The re- 
sults are summarized in Tab. ^ and plotted as 
function of Ecm (open dots) together with the 
QCD expectation in Fig. 

3. tts from Event Shapes Distributions 

From event shape distributions as is deter- 
mined by fitting an as dependent QCD prediction 
(Ofa?), pure NLLA or combined e'(a^)±NLLA 
||8|-|l0||) folded with a hadronization or a power 
correction to the data. The fitranges used for the 
different QCD predictions are shown in Fig. |^. 
The ranges for pure NLLA and 0{a1) fits are 
chosen to be distinct, so that the results are sta- 
tistically uncorrelated. 

3.1. Using Hadronization Models 

As an example of a Monte Carlo based model 
in this section the hadronization correction is cal- 
culated using the JETSET PS model (Version 7.4 
tuned to DELPHI data The QCD predic- 

tion is multiplied in each bin by the correction 
factor (the model prediction of the ratio of hadron 
over parton level in that bin). 

For 0{a1)-fits the renormalization scale /i was 
fitted to the LEPl data for both observables 
individually and then fixed for the fits to the high 
energy distributions. This takes advantage of op- 
timized scales. In contrast for NLLA and com- 
bined fits /J, was set equal to -Ecm, so that these 
results can be compared to other experiments. 

The systematic errors were obtained by fitting 



NLLA+C'(a^) 
(In i?,-schcmc) 



0.05 0.10 
NLLA 



0.35 



0.5 

1 - T 



NLLA+C'(a^) 
(In -R-schcmc) 



.03 .06 



0.25 



0.5 



NLLA 



OK) 



Ml/E^ 



Figure 1. Fit ranges chosen for fitting as from 
different QCD predictions of 1 - T and MljEl-^^ 
distributions. 



several 1— T and M"^/ E'^-^^ distributions evaluated 
by applying different cuts. The mean deviation 
from the central value is used as systematic error. 
Scale errors were taken from previous DELPHI 
publications p^ . p^ . 

The as values evaluated from the distributions 
are given in Tab. |3| and shown as a function of 
Ecxa in Fig- ||b. The results agree well with those 
measured from the event shape means. 

3.2. Using Povifer Corrections 

Dokshitzer and Webber recently extended their 
model of power corrections to event shape distri- 
butions For distributions the power correc- 
tions are applied by shifting the complete predic- 
tion: 

FU)^—^^F,,,,{f-5f) (3) 

CTtot «/ 

where the size of the shift is given by the size of 
the power correction to the corresponding mean 
{^f — (/pow))- The perturbative prediction Fpert 
can now be any of the available QCD predictions: 
0(as), NLLA or any combined scheme. 

To measure as at a single energy, again one 
first has to determine ag. This is done in a com- 
bined fit of ao and as to measurements [||j|J^ at 
many different energies. See Fig. c and Tab. ^. 
For / S^jg some of the lowest energies were re- 
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Figure 2. a) Measured mean values of (1 — T) and {M?/E^^^'j as a function of the centre of mass energy. 
The sohd lines present the results of the fits with Eqs. the dashed lines show the perturbative part. 

b) Energy dependence of as- The errors shown are statistical only. The band shows the QCD expectation 
of extrapolating the world average to other energies. 

c) as determined using Dokshitzer and Webber power corrections with 1 — T distribution. 



moved from the fit, because the shifted QCD pre- 
diction was no longer defined at the lower end of 
the fit range. 

To compare the quality of these fits to a 
hadronization model based description, a simul- 
taneous fit of as to these data using the JETSET 
PS based hadronization correction is done. The 
results are plotted in Fig. ^d, d and summarized 
in Tab. |. 

Using these ao results to determine as from 
individual high energy 1 — T distributions lead to 
systematic deviations from the QCD expectation 
(See Fig. ^). The deviation is different for NLLA 
and 0{al) (which use different fit ranges) and 
opposite in sign. 

To investigate this problem, the Dokshitzer and 
Webber power correction term was replaced by a 
simple ansatz: 



A combined fit of Ci, C2 to the complete set 
of distributions probes the need of the quadratic 
power correction term, as was fixed to 0.120 here. 
It is found, that for the fit ranges used in 0{al) 
and NLLA fits of Thrust the correction arising 
from the quadratic power term is not negligable 
(See Tab. ||). The small result of C2 for the com- 
bined result thus may be accidental. 



4. Summary/Conclusions 

Several methods of determining as at high en- 
ergies from event shapes were discussed. 

The results from means using power corrections 
are compatible with those from distributions us- 
ing standard hadronization corrections. Results 
from power corrections yield smaller scale errors. 

Power corrections to event shape distributions 
are more tricky. 
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xVndf = 416/108 



• DELPHI 172GeV 
I DELPHI 161GeV 

DELPHI 133GeV 
DELPHI 91.2GeV 
> ALEPH 91.2GeV 
o OPAL 91.2GeV 
4 SLD 91.2GeV 
» Amy SSGeV 

* Tasso 44GeV 
Tasso 34GeV 

* HRS 29GeV 

• Mark II 29GeV 
■ Tasso 22GeV 

Tasso 14GeV 




Fitrange 



» Tasso 22GeV 
• Tasso 14GeV 



0.05 0.1 0.15 0.2 



Power Correction 



JETSET 



Figure 3. 0{a'l) fit of as (and ao) to 1 — T and M^/E'^-^ distributions measured at different energies. In 
plot a,c the power ansatz Eq. ^, in b,d JETSET PS based correction is used to include non-perturbative 
contributions. The distributions including dotted lines were not included in the fit. 
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Table 3 

ag as obtained from distributions using JETSET. 



133GcV 


Q^s (-Ecm)^ cxp. ± scale 


NLLA 

CombiiKxl 


0.1154 ± 0.0072 ±0.006 
0.1196 ± 0.0129 ±0.005 

().11<S2 ± O.OOcS;-! iO.OOG 


leiGeV 


as(-Ecm)± stat. ± sys. ± scale 


NLLA 
Combined 


0.1084 ±0.0088±0.0049±0.006 
0.1056 ±0.0096±0.0043±0.005 
0.1120 ±0.0074±0.0032±0.006 


172GeV 


as(i?cm)± stat. ± sys. ± scale 


OR) 
NLLA 

Combined 


0.1151 ±0.0082±0.0058±0.005 

0.0976 ±0.0097±0.0046±0.005 
0.1097 ±0.0076±0.0029±0.005 



Table 4 

ao as determined from event shape distributions. 
The error is the statistical error from the fit. 



1 -T 


ao 


xVndf 


OR) 

NLLA 
Combined 


0.550 ±0.012 
0.475 ±0.006 
0.491 ± 0.005 


120/107 
81/26 
377/137 




ao 


X^/ndf 


OR) 

NLLA 
Combined 


0.698 ±0.014 
0.480 ±0.012 
0.554 ±0.005 


83/49 
3.2/7 
285/53 



Table 5 

Results of fitting the simple power ansatz to event 

shape distributions. o.c,=0.120. 



1 -T 


Ci[GeV] 


CalGeV^] 


O(a^) 


1.23 ±0.09 


-10.4 ±1.3 


NLLA 


0.27 ±0.02 


+10.0 ±0.6 


Combined 


0.73 ±0.02 


- 0.48 ±0.61 




Ci[GeV] 


C2[GeV^] 


0{ai) 


0.27 ±0.09 


+14.8 ±6.9 


NLLA 


not converging 


Combined 


0.70 ±0.02 


- 4.5 ±6.9 
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